ENVELOPPE SPHERIQUE SOUMISE A UNE PRESSION INTERIEURE 
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On considere une enveloppe spherique, homogene, de rayon interieur a, 
de rayon exterieur b. Le materiau qui la constitue est elastique parfaitement 
plastique, a elasticite lineaire isotrope, ayant pour critere de plasticite le 
critere de von Mises ou celui de Tresca. 

Partant de l’etat initial naturel, on soumet cette sphere a une pression 
interieure normale uniforme p que Ton fait croitre a partir de 0 (Fig.l). 



Figure 1 : Geometrie et chargement applique a la sphere sous pression 



Figure 2 : Progression de la zone plastique a partir de la surface interieure 


1 Analyse elastique 


1.1 Donner la solution (champs des contraintes et des deplacements) en 
elasticite. 

Le volume etudie est a symetrie spherique, constitue d’un materiau 


homogene et isotrope ; les conditions aux limites possedent aussi la 
symetrie spherique. On est done amene a chercher une solution du 
probleme dans un systeme de coordonnees spheriques r, 9, <]), tel que 
les champs de deplacement, de contrainte et de deformation soient 
respectivement de la forme : 

u r = h(r) uq = u§ = 0 

G rr = /i (r) dee = G^ = gi(r) G,-e = G r( j, = cj^q = 0 

t'/T — fl( r ) £qq — — Slij ) £>6 — ^r(|) — ^0 — 0 


Les equations d’equilibre se reduisent a : 
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Les conditions aux limites statiques ont la forme : 

a IT (r = a) = -p, a IT (r = b) = 0 

. . . „ , . du r 

Les equations cinematiques ont la forme : z rr = , 

dr 

La loi d’elasticite de Flooke donne : 
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Etrr = [G rr -2vGee] ££06 = [Gee(l - v) -vg, t ] 

En remplaqant les deformations par leur expression en fonction des 
deplacements, on obtient pour les contraintes les relations suivantes, 
X designant le coefficient de Lame (k = Ev/(1 — 2v ) / ( 1 + v) ) : 
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En substituant ces deux relations dans les equations d’equilibre, on 


obtient 1’ equation differentielle suivante : 
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La solution de cette equation est : u r = C\r-\ — T 

r z 

En remplaqant la valeur de u r dans les expressions precedentes, on 
obtient : 
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Les constantes C\ et C 2 s’obtiennent a partir des conditions aux 
limites : 


o rr (r = b) = 0 => C 2 = 
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Finalement, on obtient : 
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1.2 Determiner la charge limite d’elasticite P e de la sphere sous pression 
pour les criteres de von Mises et Tresca. 

Le critere de plasticite de Tresca, cornme celui de von Mises, est 
independant de la pression moyenne. On peut done l’ecrire en 
remplaqant le tenseur g par la somme de g et d’un tenseur spherique. 
Ici, si Ton ajoute a g le tenseur — Geel, on obtient un tenseur uniaxial, 
d’unique composante non nulle c rr — c*ee- D’apres les formules 
precedentes, tant que l’enveloppe spherique reste elastique, on a : 

3 a 3 b 3 

O r r ~ <^90 = yi T “ T P 

Le critere de plasticite est atteint lorsque (a,,- — Gee), fonction 
decroissante de p, devient egale a la limite d’elasticite — c y en 
compression simple. Le premier point plastique apparait done en 
r = a et lorsque la pression p atteint la valeur P e , limite d’elasticite 
initiale de la sphere sous pression : 



2 Analyse elasto-plastique 

2.1 Donner la solution (champs des contraintes et des deplacements ) en 
elasto-plasticite. Verifier que la zone plastique se developpe a partir 
de la face interne de la sphere creuse (Fig. 2). Donner la relation 
entre le rayon de la zone plastique c et la pression p. Detenniner la 
pression limite conduisant a la rupture par deformation excessive, 

Pp- 

Lorsque la pression interne p croit au-dela de la valeur P e , cornme 
le premier point plastique est apparu sur la face interieure de 
l’enveloppe, il est normal de supposer qu’une zone plastique se 
developpe a partir de cette face, et occupe un volume a < r < c, ou c 
est une fonction de p. La zone c < r < b est alors elastique. Le vecteur 


contrainte sur une facette normale a l’axe r prend la meme valeur 
dans la zone elastique et dans la zone plastique, a la traversee de la 
surface r = c. Sur cette surface la contrainte normale est alors egale 
en valeur absolue a la limite d’elasticite initiale d’une sphere creuse 
de rayon interieur c, de rayon exterieur b, soumise a une pression 
interne, soit : 

I H) 

Les contraintes dans la zone elastique sont done donnees par les 
equations precedentes dans lesquelles on remplace a par c et p par 
— o rr (c ) ; dans cette zone : 


Me) = ( 1 - 7^ ) 


O rr = ~ 


2 c 3 

3 P 


<^ee = 


2 c 3 

3P 


1 + 


- 1 a 


b 3 


2 r 3 


u r = 


2 c 3 

3 e¥ 


(1 -2v)r+(l+v) 


b^ 

2 r 


Etudions maintenant la zone plastique a < r < c. Pour y determiner 
les contraintes, on dispose des equations d’equilibre et du critere de 
plasticite, verifies en tout point, soit : 

do , r 2 \ n 

— — + -(Crr-Coo) =0 

dr r 

<5rr — <^ee = — Cty 

En combinant ces deux equations, on obtient successivement : 

do rr 2 . , 

— Oy = 0, , o rr = 2 Oy ln(r) + C 3 

dr r 

La determination de la constante d’integration C 3 s’effectue en r = c, 
en utilisant la continuite de la composante o rr : 


2o y ln(c) +C 3 = 
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Finalement, dans la zone plastique, on trouve : 
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Ces contraintes dependent du parametre c, dont il faut done 
determiner revolution en fonction de la pression p. Dans la zone 
plastique, pour r = a, on a : 


o, T (r = a) = —p 
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La transformation de la sphere etant supposee infinitesimale, a et b 
sont des constantes. La derivation de p par rapport a c donne alors : 


dp _ 2 o y c 3 

dc c ft 3 


Ce terme est toujours positif. Le rayon c de la zone plastique croit 
done constamment avec p ; ce resultat est coherent avec l’hypothcse 
que nous avons faite que la zone plastique se developpe a partir de 
la face interne de la sphere creuse. Le rayon exterieur de cette zone 
atteint la valeur b lorsque p atteint la pression limite P p : 

P p = 2a y ln 



2.2 Determiner les deformations plastiques et leurs vitesses. 

S’il est possible, a partir de la solution en contrainte obtenue 

a la question precedente de construire, en utilisant la loi de 
comportement, un champ de deplacement qui soit compatible avec 
les liaisons (cinematiquement admissible), la solution trouvee est 
unique. 

Conservant l’hypo these de symetrie spherique, on calcule le 
deplacement radial dans la zone plastique. Comme la deformation 
plastique ne produit pas de variation de volume du materiau, cette 
variation n’est due qu’a la partie elastique de la deformation, soit : 

+ 2e e e = — — [o„ + 2a ee ] 


On obtient done ainsi l’expression de la composante radiale du 
deplacement : 
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En utilisant le fait que le deplacement radial est continu a la traversee 
de la surface r = c, on peut determiner la constante d'integration : 


Ca = (1 — v) 



On obtient finalement dans la zone plastique : 
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A partir de cette expression du deplacement, on peut calculer les 
deformations totales dans la zone plastique. On obtient alors les 
deformations plastiques par difference entre les deformations totales 
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et les deformations elastiques calculees en utilisant les formules 
donnant les contraintes. Les seules composantes non nulles sont : 


2g v 


e P rr=^( l-V) (1-3 


e ee = e W = “? ((i-v)! 1 -^)) 

Comme c est une fonction croissante de p et que le trajet de 
chargement etudie est par hypothese a «p croissant», on peut 
choisir c comme parametre de chargement. Le tenseur de vitesse de 
deformation est du type compression simple : 



2.3 Que se passe-t-il si Von effectue le trajet de charge suivant : 

0 > p m (Pm A P f >) * 0 * p m ( Pm A Pe) ? 

Si L on a sounds une sphere creuse a une pression interne 

P m > P e , et qu’on la decharge jusqu’a p = 0, les contraintes 
residuelles, presentes apres cette decharge, seront egales a la 
difference entre les contraintes calculees en elasto-plasticite et la 
solution elastique correspondant a un chargement —P m . On obtient 
alors un champ de contraintes residuelles o' : 

- dans la zone plastique (a < r <c) : 
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dans la zone elastique (c < r < b) : 
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Les equations precedentes ne sont valides que s’il n ’apparait aucune 

deformation plastique pendant la decharge. Pour que la plastification 
reapparaisse en compression, il faut traverser le domaine d’elasticite, 
et retrouver un point pour lequel : a’ n . — Gg 0 = a y . Cela se produira 
effectivement a partir du moment oil la pression maximale atteinte P m 
est superieure a 2 P e . L’ etude des variations de P e et P p en fonction de 
(b/a) montre que cela n’est possible que si la pression limite est elle- 
meme superieure a 2 P e . Ceci fournit une condition geometrique sur 
la sphere. La figure 3 illustre le fait que P p depasse 2 P e si le rapport 
(a/b) est inferieur a une valeur critique x, solution de P equation 
(4/3) ( 1 —x 3 ) = 21n(x), soit : 


a/b < x ~ 0.59 


Dans le cas ou il n’y a pas plastification a la decharge, on dit que la 
structure est adaptee. Il s’agit de regime de fonctionnement sur, qui 
est utilise dans la pratique pour les recipients sous pression : ceux-ci 
subissent avant mise en fonctionnement une operation de timbrage 
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au cours de laquelle ils sont portes a une pression superieure a la 
pression de service ulterieure. 

Si au contraire il y a replastification, des deformations plastiques 
cycliques vont se produire, avec un phenomene de fatigue plastique 
du materiau, qui conduira a la ruine de la structure aux cours des 

cycles successifs 0 — > P m — > 0 — > P m — > 

La figure 4 reproduit les variations des differentes composantes du 
tenseur des contraintes a pression maximale et apres decharge, dans 
le cas oil (a/b) = 0.75. 
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a. b. 




c. d. 


FlG. 1 - (a), (b) Profils de contraintes dans l’epaisseur du tube pour differents niveaux de pression, (c) mise en evidence de l’avancee de la zone plastique, 
(d) contraintes residuelles apres decharge. 







